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1. Motivations and aims 



Theories that are asymptotically free at high energies, and have an infrared fixed point (IRFP) 
in the renormalization group (RG) flow of their couplings are said to be inside the "conformal win- 
dow". The typical running of the coupling, with its limiting value g* at small energies, is shown 
in Fig. 1 . The existence of a conformal window has been studied analytically for supersymmetric 
theories, see e.g. Ref. [1] for a review. Quantitative analytical studies of the IR regime for the 
non-supersymmetric cases are more difficult, due to a lack of adequate tools to investigate the non- 
perturbative dynamics of gauge theories. Nevertheless interesting results have appeared in recent 
years, see e.g. Refs. [2, 3, 4, 5] for a summary. Besides its field- theoretical interest, the existence 
of a conformal window has important consequences for building models of physics beyond the 
standard model (BSM) that are based on strongly-interacting theories. 
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Figure 1: Running of the coupling as a function of the energy scale for a theory with an IRFP. At low 
energies the coupling flows to a fixed-point value g*, while the high energy behaviour is the usual one 
expected for asymptotically free theories. The scale Au corresponds approximately to the energy where the 
running starts to be dictated by asymptotic freedom. The dashed curve at low energies shows the running of 
the coupling when a fermionic mass term is switched on. 



Models of dynamical electroweak symmetry breaking (DEWSB) provide elegant reaUzations 
of the mechanism responsible for the breaking of electroweak symmetry down to electromag- 
netism, by means of the vacuum expectation value of a fermion bilinear, see Ref. [6] for a re- 
view and a list of references. In the simplest versions of DEWSB, the fermion condensate is 
produced by the strong dynamics of some new asymptotically free non-abelian gauge theory. The 
new gauge group is usually called "technicolor" (TC) [7, 8]. The massless fermions coupled to this 
new gauge field are called "technifermions"; they transform under some representation of the new 
gauge group, and their left-handed components are weak doublets, so that their condensate breaks 
electroweak symmetry in the standard model (SM). The technifermion condensate also breaks the 
global chiral symmetry of the technicolor theory, in analogy with chiral symmetry breaking in 
QCD, and therefore implies the existence of massless "technipions". The number and the dynam- 
ics of technipions depend on the details of the technicolor model. Three of these technipions are 
"Higgsed" and become the longitudinal components of the W and Z bosons. As a result, the latter 
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acquire a mass proportional to the technipion decay constant, and therefore ~ v ~ 250 GeV, 
where v denotes the Higgs condensate in the SM. This constraint sets the typical scale of the tech- 
nihadronic spectrum, Atc- The rest of the technihadron spectrum depends on the specific theory 
that is selected to mediate the new strong force. Being able to compute in the strongly-interacting 
regime of the technicolor theories is the main ingredient to understand this mechanism. 

Lattice simulations are a prime tool to investigate nonperturbative physics from first principles. 
However numerical studies will play an important role in BSM studies only if we are able to identify 
the questions that are relevant in order to achieve a quantitative understanding of DEWSB. In this 
introduction we briefly review some of the problems that technicolor model building is confronted 
to, and analyze the possibility to use numerical results to make progress in this area. 

Constraints on technicolor. While technicolor provides a natural way to generate W and Z 
masses, additional interactions must be introduced in order to describe the SM flavor sector. A sce- 
nario that has been extensively studied involves "extended technicolor" gauge interactions (ETC) [9, 
10]. The ETC gauge bosons couple to both the SM fermions and the technifermions. At some en- 
ergy Metc> larger than the TC scale, the extended gauge symmetry is then broken down to the 
residual TC gauge symmetry, which remains intact. As a result, interactions between the tech- 
nifermions and the ordinary standard model fermions are generated at low-energies. We have so 
far introduced two new sectors, a TC one and an ETC one. Their connections and the energy scales 
involved are summarized in Fig. 2. Note that below the ETC scale the TC theory and the flavor 
sector of the SM are decoupled. 
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Figure 2: A schematic representation of the technicolor (TC) and extended technicolor (ETC) theories, 
and their connections with the standard model (SM). The interactions are characterized by different scales, 
indicated on the right. Bounds on the ETC scale Metc are set by SM precision measurements. 

In current simulations, only the TC sector is simulated on the lattice. At energies of the 
order of the TC scale, the ETC-mediated interactions are described by the effective lagrangian 
obtained when the heavy ETC degrees of freedom are integrated out. This lagrangian includes 
four-fermion operators, whose couplings ai^e suppressed by inverse powers of the ETC scale. These 
are dimension-six operators involving technifermions and SM fermions; denoting the SM and TC 
fermion fields by t/A(x) and Q{x) respectively, the terms that are relevant for our discussion can be 
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schematically written as: 

1 



Qix)Qix)Y{x)Yix) , (1.1) 



-^\}f{x)\if{x)Y{x)Y{x), (1.2) 

™ETC 

where for simplicity we have suppressed the flavor, color and spin indices. 

Eq. (1.1) shows that a technifermion condensate (QQ) generates a mass term for the SM 
fermions. It is worthwhile to recall here that the quark masses in the SM are defined in a given 
renormalization scheme and at a given scale /i, see e.g. the review in Ref. [6]. When the fermion 
masses are generated via an ETC theory, this dependence is reflected in the scheme and scale 
dependence of the technifermion condensate. We are going to follow the convention used in the 
PDG and always refer to the quark masses in the MS scheme. However it should be clear that this 
choice is arbitrary, and that it is possible to convert from one scheme to the other. Below the scale 
Metc the technicolor and the flavor sectors decouple, and therefore the running of the quark masses 
is simply given by the QCD logarithmic evolution, which we neglect here. In an ETC scenario, 
the SM fermion masses aie therefore set by the value of the technifermion condensate in the MS 
scheme at the scale Metc- 
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where 7 is the mass anomalous dimension of the technicolor theory in this scheme, and the second 
equality shows explicitly the evolution of the chiral condensate. 

On the other hand, ETC interactions in Eq. (1.2) include flavor changing neutral currents 
(FCNC), which would result in sizeable contributions to K and D mixing. As noted long ago [9], 
the experimental limits on such contributions imply a lower bound on the ETC scale. For instance 
the ETC scale associated with the generation of the strange mass should be larger than 10^ TeV, 
unless some fine-tuned mechanism is in place. As a consequence of this large value of METC^ 
heavy quark masses can be generated only if the techniquark condensate is enhanced with respect 
to the naive value obtained by rescaling QCD. A detailed discussion of the fermion mass generation 
in ETC models can be found in Ref. [11]. Recent summaries on technicolor and the constraints 
imposed by current measurements can be found in Refs. [2, 12, 13]. 

It has been suggested that "walking" theories, i.e. theories where the RG evolution is very 
slow, could produce the required enhancement of the condensate [14, 15, 16, 17, 18, 19, 20]. This 
is illustrated in Eq. (1.3): if the couplings ai^e approximately constant as the energy scale varies, 
the integral of the anomalous dimension yields a power enhancement: 



^^e)U.. = (^)'(e0L,,. (1.4) 



For comparison note that, in a QCD-like theory, the condensate at the higher scale is of the order of 
(22) 1^^^ up to logarithmic connections. For a sufficiently lai^ge value of 7, the factor in Eq. (1.4) 
can yield large quai^k masses and hence ease the tension between technicolor and the flavor sector 
of the SM. Recent analyses suggest that 7 > 1 is needed to avoid a fine-tuned ETC sector [21], 
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while unitarity constraints imply that 7 < 2 [22]. For a strongly-coupled theory, the value of / can 
only be computed by nonperturbative methods. 

Furthermore "walking" theories are expected to ease the tension between TC and the con- 
straints from electroweak precision measurements [23, 24, 25]. 

A lucid list of questions that need to be addressed in order to build a phenomenologically 
successful model of DEWSB is presented in Ref. [26]. 

It is important to recall that the rate at which the coupling constants evolve, and the corre- 
sponding power enhancement, depend on the renormalization scheme, and therefore the definition 
of "walking" needs to be qualified better. The existence of an IRFP, i.e. a zero of the beta functions 
for the couplings at small energy scales, is independent of the scheme. It describes the physical 
properties of theories that are scale-invariant at large distances, where the field correlators have 
power-like behaviours characterized by the anomalous dimensions of the fields ' . A concise de- 
scription of the scheme dependent features of IRFPs is presented in Ref. [27]. 

The conformal window. In the absence of fermionic matter fields, SU(N) gauge theories are 
asymptotically free. The running of the gauge coupling g is encoded in the beta function, which 
can be computed in perturbation theory close to the Gaussian fixed point g = 0: 

I^J^S = P{8) = -Pog'-Pig' + ■■■■ (1.5) 

At one-loop in perturbation theory the effect of the fermion fields can be read from the first coeffi- 
cient 



11 , , 4 
-C2{A)--Tj,nf 



(1.6) 



(471) 

where C2{R) and Tr denote respectively the quadratic Casimir and the normalization of the gener- 
ators in the representation R. The coefficient C2(A) originates from the gluons being in the adjoint 
representation of the gauge group, and n f counts the number of Dirac fermions in the theory. At 
one loop the dependence on the fermionic representation is entirely encoded in the factor Tr. In 
discussing perturbative results, we use here the same conventions introduced in Ref. [28]. As the 
number of fermion fields is increased, jSq changes sign and asymptotic freedom is lost. We shall 
denote by «/^up the number of fermions above which asymptotic freedom is lost, cleai^ly such a 
number depends on the gauge group and the fermionic representation. 

At higher orders in perturbation theory the fermionic contribution to the running of the cou- 
pling has the potential to generate a non-trivial zero of the beta function before asymptotic freedom 
is lost. This is signalled in perturbation theory by a change of sign of the coefficient j8i . We shall 
denote by n fjo the number of fermions above which the theory exhibits an IRFR In this case the 
theory becomes scale-invariant at large distances, while the short-distance behaviour is still the one 
dictated by asymptotic freedom. As a consequence of scale invariance at large distances, the theory 
cannot be in a confining phase and chiral symmetry remains unbroken. The long-distance dynam- 
ics is governed by the IRFP's critical exponents, which determine the scaling laws in the vicinity of 
the fixed point. The range of values Uf jo <nf < nj up is known as the "conformal window". The 



'Anomalous dimensions at tiie fixed point are closely related to the critical exponents that are introduced to study 
statistical systems near criticality. 
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Banks-Zaks theories [29, 30], where Nc and nf aie arranged such that the critical coupling g* ^ 1, 
provide one working example of theories within the conformal window that can be analysed in 
perturbation theory. 

Theories near the edge of the conformal window (e.g. «/ < «/.io) have been recently put for- 
ward as candidate walking theories [31]. Alternative scenarios start from a theory inside the con- 
formal window, and deform the theory away from conformality by perturbing it with operators that 
are relevant in the IR regime [32], see also Refs. [33, 34]. In both cases, the starting point is being 
able to identify the conformal window, i.e. to be able to identify the existence of a fixed point, and 
to compute the critical exponents that characterize the relevant directions of the RG flow. We shall 
concentrate on this specific problem in this review. 

Informations about the existence of IRFPs in nonsupersymmetric gauge theories can be ob- 
tained by analytical methods like e.g. the solution of the truncated Schwinger-Dyson equation [35, 
36, 37, 38], or the usage of conjectures about the nonperturbative behaviour of quantum field the- 
ories [39]. This information is valuable to guide the numerical investigations, but robust results 
ultimately need investigations that rely on first principles; besides the lattice studies discussed be- 
low, interesting analytical results for the anomalous dimensions in a CFT have been obtained in 
Ref. [40] from first principles. 

Lattice evidence. Having set the task to identify the existence of IR fixed points, we need to 
spell out clearly what are the observables that can be evaluated by Monte Carlo methods, and what 
are the systematic errors that need to be kept under control in order to draw meaningful conclusions 
from the lattice data. 

The tools that have been developed for the numerical studies of QCD can be adapted to this 
new class of theories, yielding new and interesting informations. These tools are reviewed in Sect. 2 
below. 

In order to perform a lattice study, a particular theory has to be selected and simulated. Analyt- 
ical results can help in selecting candidate theories inside the conformal window. However the lack 
of robust analytical results means that there is no guarantee of being inside (or near) the conformal 
window before simulations are performed. As a consequence, numerical tools developed for these 
studies need to provide: 

• flexibility; so that the same codebase can be used to simulate more than just one theory. 

• efficiency; so that simulations are reasonably fast. At this stage, it is important to find the 
right balance between optimization and flexibility. 

• a wide range of observables; so that robust conclusions can be drawn, based on more than 
one observation. 

• detailed numerical benchmarks, so that the algorithmic issues can be clearly identified, and 
distinguished from the physically meaningful results. In my opinion these algorithmic stud- 
ies have not yet reached the maturity of their QCD counterparts. 

A few candidate theories have been studied recently, chosen according to the conjectured 
boundaries of the conformal window for SU(N) gauge theories with fermions in the fundamental. 
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adjoint or two-index representations. These boundaries liave been summarized in Ref. [2], whence 
the plot in Fig. 3 is extracted. As shown in this plot, theories with fermions in higher dimensional 
representations of the gauge group are expected to develop a fixed point for relatively small num- 
bers of flavors. Recent studies have focused on theories that are close to the lower edge of the 
conformal window: SU(3) with nf = 16, 12, 10,9,8,6 flavors in the fundamental representation, 
SU(2) with rif = 6 flavors in the fundamental, SU(2) with nf = 2 flavors in the adjoint represen- 
tation, and SU(3) with nf = 2 flavors in the two-index symmetric {sextet) representation. At these 
early stages of the nonperturbative studies of the conformal window it is important to try to identify 
a paradigm to guide the numerical investigations, rather than trying to get exhaustive results on one 
specific theory. 



SU(N) Phase Diagram 
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Figure 3: Boundaries of the conformal window for SU(N) gauge theories with species of Dirac fermions. 
The four bands represent respectively fermions in the fundamental (Fund), adjoint (A) and two-index sym- 
metric and antisymmetric (2S,2A) representations. The upper limit of each band corresponds to the number 
of flavors where asymptotic freedom is lost, as obtained from one-loop perturbative computations. The 
lower limit of each band yields the number of flavors above which the theories develop an IR fixed point. 
The location of these lower limits relies upon assumptions about the nonperturbative dynamics of the theo- 
ries. Lattice simulations can provide first-principle evidence in favour (or against) this picture, and compute 
the critical exponents that characterize the fixed points. Figure courtesy of F. Sannino. 

2. Tools 

Numerical tools that were originally designed for investigating lattice QCD have been used in 
order to identify the existence of IRFPs. We describe briefly the main ideas, the observables that 
ai^e used in the different approaches, and their expected behaviour in the presence of an IRFR For 
each case we try to emphasize the sources of systematic en^ors that need to be kept under control 
in order to draw robust conclusions from numerical data. 

2.1 Phase structure of the lattice theories. 

Lattice simulations are performed by discretizing the action of a given theory on a Euclidean 
space-time lattice. At weak coupling the RG flow can be computed perturbatively, and the relevant 
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parameters are easily identified. For an asymptotically-free gauge theory, ^ = is an UV fixed 
point that defines the usual continuum limit of the lattice theory. The IRFP that we are seeking is 
a fixed point on the massless renormalized trajectory that originates from the continuum limit. As 
the bare coupling is increased, the lattice theory may have a complicated phase structure, with new 
fixed points appearing, which are not necessarily related to the continuum limit. 

In particular there can be bulk phase transitions in the lattice theory at strong coupling. These 
are lattice artefacts and can obscure the continuum physics that we are interested in. Mapping the 
phase structure of the lattice theory is important in order to be able to identify the genuine features 
of the continuum limit. 

Finite volume effects also play an important role in determining the long-distance behaviour of 
a gauge theory. The inverse size of the lattice 1/L is a relevant coupling in the IR, which introduces 
an explicit length scale: we can only probe energies above 1 /L when performing lattice simulations 
in a finite box, while the largest scaling factor that we can probe in a simulation is ^max — L,/ CI. 

There are two potential effects to take into account. On the one hand, as the box size is 
decreased, the IR cutoff can potentially change the behaviour of the gauge theory and drive the 
theory into different phases, characterized by the breaking of centre symmetries [41, 42]. It is 
important to study the impact of these transitions on the measured observables. On the other hand, 
even when the box is large enough and the system is in the correct phase, 1/L determines the 
size of the scaling violations. These can be readily included in scaling laws determined from RG 
analyses [43, 44, 45, 46, 47]. 

2.2 Spectral studies 

Informations about the mass spectrum and the decay constants of a given gauge theory, can 
be obtained from the numerical analysis of two-point correlators, using the standard techniques 
developed for QCD. 

Simulations are performed at a finite value of the mass, and in a finite volume. Both these 
quantities break conformal invariance at large distances. The signatures of an IRFP in this frame- 
work are the scaling laws of physical observables as functions of the mass and volume, as it is often 
the case when studying critical phenomena. Using the scaling behaviour of lattice observables to 
identify IRFP was proposed in Refs. [48, 49, 43, 50]. Scaling laws are derived using the RG equa- 
tions for field correlators [43, 50]. A detailed analysis of the scaling of lattice observables was 
discussed in Refs. [46, 47]. Scaling laws in the vicinity of a fixed point were already considered in 
searches for UV fixed points in lower-dimensional field theories. For a summary of results in the 
case of three-dimensional four-fermi theories, see e.g. Ref. [51] and references therein. 

The masses of the states in the technihadron spectrum scale as: 

MH~mi/>'"', (2.1) 

where = 1 + 7* is the critical exponent associated to the fermion mass, and is the (scheme- 
independent) value of the anomalous dimension at the fixed point. 

A generic operator O with appropriate quantum numbers couples to the states of the theory 
with strength denoted Ge for a scalar state, and Fy for a vector one. The scaling of these cou- 
plings, which are related to the decay constant of the state H, can also be derived from the scaling 
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behaviour of two-point correlators. Following the notation introduced in Ref. [46] the scaling be- 
haviour of the decay constants can be written as: 

G ~ m^'o , F r^m'^" , (2.2) 

with the exponents summarized in Tab. 2.2. Note that the scaling relations hold for any state in 
the spectrum. In particular for the pseudoscalar meson states, they yield a modified Banks-Casher 
relation. Scaling laws in the vicinity of a fixed point have also been discussed in the context of the 
functional RG, see e.g. Ref. [52, 53]. 
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Table 1: Scaling laws, G[F] ^ m''°l^l for decay constants. The symbol y„, = \ + denotes the scaling di- 
mension of the mass and h.o = do + Jo is the scaling dimension of the operator O. The symbol a denotes the 
adjoint flavour index, and A" are the generators normalized as tr[A"A''] = 25"''. No such simple expression 
exists for the axial singlet current because of the chiral anomaly [46] . 

It is important to bear in mind that scaling laws are asymptotic formulae that are obtained by 
hneaiizing the RG equations in a neighbourhood of the fixed point at m = and 1/L = 0. As the 
system moves towards larger values of m and 1/L, the corrections to scaling become important 
and eventually obscure the power-law scaling dictated by the IRFR The critical exponents can be 
extracted from the power-law scaling of the spectrum only from simulations at small mass. 

The distinctive feature in the spectrum of an IR-conformal theory is the lack of spontaneous 
chiral symmetry breaking. As a consequence, there are no Goldstone bosons in the theory. As 
the fermion mass is sent to zero all states in the spectrum become massless, see e.g. Ref. [54] 
for a realization of this scenario. This is at odds with the chirally-broken scenario, where there 
is a pai^ametric separation between the pseudoscalar Goldstone bosons and the rest of the massive 
spectrum. Any effective theory that describes the low-energy dynamics of an IR-conformal system 
must take all the light degrees of freedom into account. 

As simulations move to smaller fermion masses, the physical size of the lattice must be in- 
creased in order to avoid finite-size effects. In particular a chirally broken theory at small mass 
and fixed physical volume can enter the so-called 5-regime, where the spectrum is determined by 
the rotator states of the chiral condensate [55]. This regime corresponds to Lt » Ly, FLs > 1, and 
MLg <^ 1 ; the spectrum of the lowest states no longer scales as expected for the Goldstone bosons of 
a chirally broken theory. In this case deviations from the GMOR scaling are not a signal of restora- 
tion of chiral symmetry. Similarly if ~ the system will be driven in the e regime, which can 
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be seen as the high-temperature limit of the 5-regime. Great care must be exercised in the inter- 
pretation of lattice data, as these finite size effects can easily be mistaken for signals of conformal 
behaviour. This has been clearly pointed out in Ref. [56]. The 5 regime can be distinguished from 
a conformal one if the power-law scaling corresponding to the latter is observed. 

Last but not least, it is important to verify that the simulated masses do not correspond to the 
"heavy quark" limit of a chirally broken theory. This can be achieved by comparing the mesonic 
and gluonic sectors of the spectrum as discussed in the next section, where we report the numerical 
results. 

2.3 SF studies 

The nonperturbative running of the gauge coupling and fermion mass can be studied using the 
Schrodinger functional scheme [57, 58]. The running coupling at the scale 1/L is defined on 
a hypercubic lattice of size L, with boundary conditions chosen to impose a background chromo- 
electric field, which depends on a parameter rj . The renormalized coupling is defined as a measure 
of the response of the system to changes in the background chromoelectric field: 



where S is the action of the Schrodinger functional, and the constant k is chosen such that = g^ to 
leading order in perturbation theory. Eq. (2.3) defines a nonperturbative coupling, which depends 
on only one scale, the size of the system L, and can be evaluated numerically. 

To measure the running of the quai^k mass, we calculate the pseudoscalar density renormalisa- 
tion constant Zp. Following Ref. [59], Zp is defined by: 

Zp{L) = y^i/fp{L/2) , (2.4) 

where f\ and fp are the correlation functions involving the boundary fermion fields ^ and 

/i = -l/UL" ! d\d\d^yd\{t {u)YsT%'{v)X{y)Y,x%{z)) , (2.5) 

/p(xo) = -1/12 j d'yd\.{W{xo)r5r''v{xoKiy)r5T%{z)) . (2.6) 

For the details of the Schrodinger functional setup we refer the reader to the original publications. 

The running of the coupling as the scale is varied by a factor s is encoded in the step scaling 
function a{u,s) as 

J:{u,s,a/L) = f{go,sL/a)\_2(^g^j^^^^^^ , (2.7) 
g{u,s) = lim L{u,s,a/L) , (2.8) 

a/L^O 

as described in Ref. [58]. The function a(M, 5) is the continuum extrapolation of £(m, 5, a/L) which 
is obtained from numerical simulations at various a/L values and fixed u. The step scaling function 
encodes the same information as the j8 function. The relation between the two functions for a 
generic rescaling of lengths by a factor s is given by: 

-2\ogs= (2.9) 

Ju ^Jxp{y/X) 
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The step scaling function can be computed at a given order in perturbation theory by using the 
analytic expression for the perturbative j8 function, and solving Eq. (2.9) for a{u,s). It can be seen 
from the definition of a{u,s) in Eq. (2.8) that an IRFP con^esponds to a{u,s) = u. 
The lattice step scaling function for the mass is defined as: 



Lp[u,s,a/L) — 



Zp{go,L/a) 

the mass step scaling function in the continuum limit, ap{u,s), is given by 



; (2.10) 

f{L)=u 



ap{u,s) = \imLp{u,s,a/L) . (2.11) 



The mass step scaling function is related to the mass anomalous dimension (see e.g. Ref. [60]): 

(2.12) 



„ , (4/(2ft)) 



^^^/^rW do 



Op[u = exp 

o{u)J yj^, \p{x) Pox 

In the vicinity of an IRFP the relation between Op and 7 simplifies: 

™('^/*) dm fi^l' dq 

— = -7* — , (2.13) 
mill) m J pi q 

and hence: 

\og\op{s,u)\ = -yjogi. (2.14) 
We can therefore define an estimator 

^(„)^_12IM«)1. (2.15, 

which yields the value of the anomalous dimension at the fixed point. Away from the fixed point 7 
will deviate from the anomalous dimension, with the discrepancy becoming larger as the anomalous 
dimension develops a sizeable dependence on the energy scale. 

The step scaling functions need to be extrapolated to the continuum limit in order to disentan- 
gle the running of the couplings from the lattice artefacts that affect the measured observables. This 
extrapolation is very delicate for a{u,s). In the vicinity of a fixed point the running is by definition 
very slow. Therefore a very high accuracy in the extrapolation is needed in order to resolve the 
physically meaningful signal. The systematic uncertainties make it difficult to locate the value of 
the critical coupling with sufficient accuracy. In the next section, we shall discuss in more detail 
the propagation of errors and their impact on the physically interesting results. 

2.4 Potential schemes 

Finally the running of the coupling can be studied by defining a nonperturbative coupling from 
the potential between static charges computed numerically [61]. 
The coupling is defined as: 

gl{Lo,R/Lo,a/Lo) = ^{R/afx{R/a,Lo/a) , (2.16) 
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where ;^ is a Creutz ratio, Lq is the size of the lattice, and R is the size of the Wilson loops used to 
construct the Creutz ratios. The size of the lattice Lq sets the scale at which the coupling is defined, 
while the ratio R/Lq defines the renormalization scheme. 

A step scaling function can be defined in close analogy to the one defined for the Schrodinger 
functional: 

i:„{u,R/LQ,a/Lo) = gl{bLQ,R/LQ,a/Lo) , (2.17) 
u = gl{Lo,R/Lo,a/Lo), (2.18) 

o„(u) = lim'Lw(u,R/LQ,a/Lo) . (2.19) 

Once again the extrapolation to the continuum limit in Eq. (2. 19) is necessary in order to avoid 
contaminations from lattice artefacts. 

Ref. [62] presents a variation on the same theme. The potential between static charges is 
measured by lattice simulations, and is compared with the results obtained from integrating the 
running coupling: 

V{R) -V{Ro) = QiR) £dR' ^ , (2.20) 

where the perturbative expansion is used as an input for a (R) in the integral. This method allows 
one to compare the nonperturbative running with the perturbative expectation. 

2.5 MCRG 

Monte Carlo RenormaUzation Group (MCRG) methods were developed in order to study the 
coupling flow in both spin and gauge models. In particular, the 2-lattice matching has proved to be 
useful in pure Yang -Mills theories [63, 64, 65]. 

The basic idea is to follow the RG flow of the bare couplings under blocking transformations 
that integrate out the UV degrees of freedom. With each blocking step, changing the scale by a 
factor s, the flow drives the couplings towards a lower-dimensional manifold, whose dimension is 
given by the number of relevant couplings. The distance from this manifold goes as a power of s. 

2.5.1 Two-lattice matching procedure 

Let us consider first for simplicity a theory that has only one relevant parameter flowing out 
from an UV fixed point; this is the common situation in pure gauge theories. In this case the 
RG trajectories converge towards a one-dimensional renormalized trajectory (RT). Starting from 
a value g of the bare coupling at the cutoff scale, after n steps the Wilsonian action is described 
by some point in pai^ameter space, for a sufficiently large n this point is close to the RT. This is 
represented by the circles and the blue curve in Fig. 4, for ?i = 4 steps of blocking starting from the 
value g = K. This point is matched by the point g' = K' such that the RG flow from K' ends up at the 
same point as the previous flow after only « — 1 steps. This flow is represented by the diamonds on 
the red trajectory for n = 3. Since the two trajectories end at the same point, the lattice coiTclation 
length at the endpoint of the two flows must be the same for both theories: | {g)/s" = | {g')/s"^^, 
and thus: 

a{g')=sa{g). (2.21) 
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To identify such a pair of couplings, we need to show that after n and {n — I) RG steps re- 
spectively their actions are identical. Explicitly calculating the actions would be complicated, but 
instead the gauge configurations themselves can be blocked; showing that the expectation values of 
all observables on these gauge configurations agree is equivalent to directly comparing the actions 
that generated them. 

This procedure identifies, for each n, a pair of bare gauge couplings {g,g'), or equivalently 
(j8,j3') where j8 = 2N/g^, with lattice correlation lengths that differ by a factor s, |' = ^/s. In 
the limit n — oo, it is customary to define the step scaling function for the bare coupling: Aj8 = 
^ — ^' = Sh{^',s). This is the analogue of the Schrodinger Functional step scaling function for 
the renormalised coupling, o{u,s), and in the UV limit where ^ g}) = 2N / j8, there is a simple 
relation between the two: 

^^ = ^-1. (2.22) 
p u 

Clearly an IRFP is found when si,{^*;s) = 0, while Sj, is expected to remain positive for a QCD-like 
theory. 

There is a degree of arbitrariness in choosing the blocking transformation, see e.g. Ref. [66]. 
Recent studies have used: 



V„,p = Proj 



(2.23) 



where a is a free parameter, which can be varied to optimise the transformation. Another possible 
choice for blocking is to perform a so-called HYP-smearing [67, 68]. Changing the blocking 
ti'ansformation changes the location of the fixed point, and the rate of convergence towards the RT. 
Ideally it should be chosen such that: (i) All observables yield the same {g,g') pairs for a given 
number of blocking steps n. Deviations are a measure of the systematic eiTor from not being at 
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exactly the same point along the RT; (ii) consecutive blocking steps predict the same matching 
coupling, i.e. for a given g, the matching coupling g' should be the same for all n. Deviations show 
that the irrelevant couplings still have sizeable effects. MCRG yields robust information only if 
these systematic errors are under control. 



3. Results 2010 



We shall summarize here the latest results at the time of the Lattice Conference. Previous 
studies have been summarized in the plenary talks at the Lattice Conferences in 2008 and 2009 [69, 
70]. 

3.1 SU(3) with fundamental fermions 

Starting from the upper end of the conformal window «/ = 16, several theories have been 
studied with decreasing numbers of fermions. 

At Uf = 16 the two-loop beta function predicts an IRFP at weak coupling, g^ w 0.5. MCRG 
studies of this theory [67] have found a negative step scaling function Sh, in agreement with the 
hypothesis that the theory is indeed inside the conformal window. The bare step scaling function is 
shown to the left of Fig. 5. Different types of blocking yield different locations of the zero of the 
step scaling function, as expected since the position of the fixed point depends on the renormaliza- 
tion (or blocking) scheme. 



o: ORIG 
X: HYP 



0.4 



N,= 12 



o: ORIG 
X: HYP 
i:;: HYP2 



T X* * f I 



Figure 5: Bare step scaling function for «y = 16 (left) and «^ = 12 (right) using the ORIG (o), and the HYP 
(x) blocking procedures; details about the blocking procedure can be found in Ref. [67]. 



At Hf = \2 the situation is less clear. MCRG studies suggest a slow running of the coupling, 
but do not provide a clear cut answer about the existence of a fixed point. The bare step scaling is 
shown on the right of Fig. 4. The main limitation is the fact that these simulations require to run the 
code at exceedingly large values of the bai^e coupling [67]. Larger lattices could help to improve 
these results by allowing more blocking steps. 

Studies of the running of the coupling defined through potential schemes do not show any sign 
of a conformal fixed point in this case. We refer the reader to the talks of Itou and Holland in these 
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Proceedings [71, 72] for more details on these studies, which disagree with the results presented in 
Refs. [73, 74]. 

Spectral studies for this theory have not reached consensus yet. Some are consistent with the 
spectrum of a confining theory [75, 76, 56]. A summary of these results is reported in Fig. 6. 
However other studies fit well the hypothesis that the theory is in a chirally symmetric regime [77, 
78, 79, 80]. As discussed in the previous section, there ai^e a number of systematic errors that ai^e 
likely to obscure the physically meaningful results, and more extensive simulations will be required 
to settle this issue. 



Lattice Higgs Collaboration 



New Nf=12 NLO beta=2.2 chiral analysis in p-regime: 




0.01 0.02 0.03 0.04 



Consistent with chiral symmetry breaking 



Figure 6: Spectrum studies for SU(3) with «/ = 12. The upper plots show the pseudoscalar masses. The 
lower plots display, from left to right, the GMOR relation, the pion decay constant, and the chiral condensate. 
Lines correspond to fits to chiral perturbation theory [75]. 

As is further decreased, results for «j = 10,9,8 indicate these theories are already below 
the edge of the conformal window [67, 56]. 

Recent results, which appeared after the Lattice conference, suggest instead that the theory 
with «/ = 10 lies inside the conformal window [81]. 

The LSD collaboration has been studying the «/ = 6 theory. This theory is expected to be in 
the confining regime, sufficiently close to the edge of the conformal window to display a walking 
behaviour. An enhancement of the chiral condensate by a factor of 2 at the cutoff scale, and 
a possible trend towards parity doubling for the mesonic states have been identified [82]. The 
phenomenological imphcations of these findings need to be studied in more detail. 

3.2 SU(2) with adjoint fermions 

The SU(2) gauge theory with fif = 2 flavors in the adjoint representation has been investigated 
by several groups using different methods. In this case all simulations seem to indicate that the 
theory is inside the conformal window. 
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Recent studies are reported in Refs. [83, 28, 84, 85, 86, 87, 88, 50, 89]. Simulations so far 
have been performed with non-improved Wilson fermions; the phase diagram for the lattice theory 
has been mapped carefully in Refs. [87, 85], where a bulk phase transition was found and the 
region connected to continuum physics has been identified. Simulations have been performed 
trying to reach the small mass regime while preserving the hierarchy of scales required to control 
the systematic errors: 

a 



<^ amps <■ 



< 1. 



(3.1) 



where L is the lattice size, mps is the lightest mass in the mesonic spectrum, and ro is the Sommer 
radius. Satisfying the above inequalities at small fermion masses becomes very rapidly a CPU- 
intensive task. New results for the spectrum were presented at this Conference [90, 91, 92]. Fig. 7 
summarizes the most striking features observed in the spectrum, namely the near-degeneracy of 
the pseudoscalar and the vector meson mass, which persists at the smallest masses explored so far, 
and the large ratio of the pseudoscalar mass to the string tension. Both behaviours are at odds with 
the expected behaviour in a theory where chiral symmetry is spontaneously broken. In the latter 
case, the pseudoscalar mass goes to zero, while the other two quantities remain finite, thus yielding 
respectively a divergent and a vanishing ratio for the quantities in Fig. 7. Another interesting aspect 
of the spectrum study presented in Ref. [93] is the hierarchy between the glueball and the mesonic 
states, with the former being lighter than the latter. The mass of the glueballs scales with the 
fermion mass m, indicating that the system is not in the heavy mass regime. This result suggests 
that the light glueball states need to be included in any effective lagrangian description of TC low- 
energy dynamics. Note that these results have been obtained at a single value of the lattice spacing 
and therefore the size of lattice artefacts cannot be estimated properly. 
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Figure 7: Vector meson to pseudoscalar meson mass ratio (left); notice that the two states remain degenerate 
when the mass is decreased. Ratio of the pseudoscalar mass to the square root of the string tension (right). 
The pseusoscalar mass and the string tension vanish at the same rate, yielding a finite ratio. 



Finite size scaling can be tested by rescaling data obtained on lattices of varying size. Assum- 
ing the existence of an IRFP, data from different lattices should fall on a universal curve, e.g. for 
the pseudoscalar decay constant: 

LFps = ^{V'"m). (3.2) 
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Finite size scaling curves are displayed in Fig. 8, where different values of the scaling exponent 
ym = 1 + 7* are used for rescaling the data. The plot indicates that the numerical results are consis- 
tent with the existence of an IRFP with a low value for y,„. 
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Figure 8: Finite size scaling curves for the pseudoscalar decay constant. The four plots, clockwise from the 
top left corner, correspond to scaling exponents y* = 0.1,0.2,0.5,0.4 respectively. 

Schrodinger functional studies for this theory show some evidence in favour of the existence 
of a conformal fixed point. The running of the gauge coupling was studied in Refs. [94, 27]. The 
results are summarized in Fig. 9. The plot on the left shows the continuum step scaling function, 
an IRFP is signalled by the condition o{u)/u = 1. The plot on the right is simply the difference 
of the renormalized couplings at scales L\ and L2; with the latter choice of variables the IRFP is 
identified by the difference of the two couplings, A(Li ,L2), changing sign. Note that in this second 
way of presenting the data, the continuum limit is not taken, but the lattices with larger values of 
L\ are closer to the continuum limit. In both cases, we see that the data are compatible with the 
existence of a tixed point. However the current error on the data (especially when performing the 
continuum extrapolation) is too large to locate precisely the critical value g*. This is not surprising 
if the theory is really conformal, or near the edge of the conformal window. When the running 
of the coupling becomes very slow, the SF simulations have to resolve a very small signal that is 
easily obscured by the statistical and systematic errors. This is a common problem of all the studies 
of the running coupling for conformal theories. 

Finally, as discussed in the previous Section, the SF can be used to study the running of 
the mass and therefore deduce the anomalous dimension 7 [27]. Data for Op{u) are displayed in 
Fig. 10. Note that the step scaling function Op is not small for a conformal theory, and therefore 
can be measured with better accuracy. Interestingly its value does not display any statistically 
significant deviation from the one-loop prediction. 

The anomalous dimension 7* can be read from the curve in Fig. 10, as long as the value of the 
gauge coupling at the fixed point is precisely known. Unfortunately, as discussed above, the latter 
cannot be located precisely; the uncertainty on the value of g* is the main source of uncertainty in 
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Figure 9: Running of the coupling constant in the SF scheme. Data are taken from Ref. [27] for the left 
plot, and Ref. [94] for the right one. The large error bands in the plot on the left are due to the continuum 
extrapolation of the step scaling function. 
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Figure 10: Running of the mass computed from the SF [27]. 



the determination of the mass anomalous dimension. With the current data the best estimate is: 

0.05 < 7* < 0.56. (3.3) 

Preliminary results using the MCRG technique seem to confirm this estimate [95]. 

Recent results confirming the picture above have recently appeared [96]. 

Results obtained with different techniques are consistent with the existence of an IRFP for 
these theory, with a small value of the anomalous dimension. The phenomenological consequences 
of such a small value need to be investigated carefully. 

3.3 SU(3) with sextet fermions 

Results have also been obtained for the SU(3) gauge theory with «y = 2 in the two-index 
symmetric representation, i.e. the sextet representation of SU(3). This model has also been pro- 
posed for phenomenological applications under the name of Next to IMinimal Walking Technicolor 
(NIVIWT). 

Results for the spectrum and the eigenvalues of the Dirac operator have been presented in 
Refs. [97, 98, 99, 100, 101, 102, 44]. 
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Results for the SF lattice step scaling function £(11) have been presented at this Confer- 
ence [103], showing evidence of a slow running of the coupling. Different discretizations yield 
statistically inconsistent estimates for the location of the fixed point coupling g*, emphasizing the 
importance of controlling the size of lattice artefacts in these studies. The difference 



Biu,2) = -l/gHL) + l/gH2L) 



(3.4) 



is displayed on the left-hand side in Fig. 1 1 , computed using thin links (o) and fat links (□) respec- 
tively. The difference between the diamond curve and the square curve shows the effect of lattice 
artefacts. 

Once again, results for the mass anomalous dimension have a much smaller relative error, as 
seen in the plot on the right of Fig. 1 1 . The estimates obtained from the SF, from finite size scaling, 
and from the eigenvalues of the Dirac operator are in broad agreement, and suggest the bound 
Y* < 0.6. 

Studies at finite temperature on the other hand indicate that the theory is outside the conformal 
window and "slow walking" [104, 105]. 
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Figure 11: Running of the coupling (left) and the mass (right) for the SU(3) theory with two flavors in the 
two-index symmetric representation [103]. 

A detailed computation of the spectrum has been performed with staggered fermions for the 
same theory, and presented at this Conference [75]. The numerical data in this case are consistent 
with chiral symmetry breaking. 

Further work is needed to reach consensus on the long-distance dynamics in this case. 



4. Outlook 

There have been numerous studies of nonperturbative dynamics beyond QCD in the last two 
years. Several theories that are expected to be at the edge of the conformal window have been 
investigated, using the tools described in Section 2. This is an interesting problem in field theory, 
and there are several intriguing data, suggesting that first evidence for the existence of fixed points 
has been found. However robust results can be obtained only by keeping the systematic eiTors 
under control. 
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B parameter for Nf = 2 sextet representation M'^ for Nf = 2 sextet representation 
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Figure 12: Chiral extrapolation of the chiral condensate and other spectral quantities. Data are consistent 
with rooted staggered chiral perturbation theory [75]. 



In particular simulations must be performed at small masses, on large volumes, trying to mini- 
mize lattice artefacts that can obscure the physically interesting features. Preliminary results about 
simulations with improved actions for the SU(2) theory have been presented in Refs. [106, 107]. At 
the same time it is worthwhile to keep looking for better observables that can yield unambiguous 
signals. 

It is important to develop a strong link with the phenomenological work and the data analyses 
at the LHC. Lattice data will have an impact on phenomenology only if the relevant questions are 
identified and answered in a quantitative way. A wishlist of interesting issues has been discussed 
by Chivukula at this Conference [108]. 
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